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ABSTRACT: Radiative antenna techniques, e.g., dipole and monopole, have been pro- 
posed for radiofrequency (RF) coil array designs in ultrahigh field MRI to obtain stronger B, 
field and higher signal-to-noise ratio (SNR) gain in the areas deep inside human head or 
body. It is known that element decoupling performance is crucial to SNR and parallel imag- 
ing ability of array coil and has been a challenging issue in radiative antenna array designs 
for MR imaging. Magnetic wall or induced current elimination (ICE) technique has proven to 
be a simple and effective way of achieving sufficient decoupling for radiative array coils 
experimentally. In this study, this decoupling technique for radiative coil array was analyzed 
theoretically and verified by a simulation study. The decoupling conditions were derived 
and obtained from the theory. By applying the predicated decoupling conditions, the isola- 
tion of two radiative elements could be improved from about —8 dB to better than —35 dB. 
The decoupling performance has also been validated by current distribution along the radia- 
tive elements and magnetic field profiles in a water phantom. © 2015 Wiley Periodicals, 
Inc. | Concepts Magn Reson Part B (Magn Reson Engineering) 45B: 183-190, 2015 
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I. INTRODUCTION to-noise ratio (SNR), higher spatial resolution, and 


better image contrast (7—5). The high static magnetic 


Ultrahigh field (i.e., 7 T and higher) magnetic reso- 
nance imaging (MRI), has demonstrated remarking 
benefits for human imaging, including higher signal- 
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fields result in the increase in the Larmor frequency 
which is usually in radiofrequency (RF) range. At 
the high frequency of 300 MHz (Larmor frequency 
of proton at 7 T) and beyond, the wavelength is sig- 
nificantly shortened, especially in high dielectric bio- 
logical samples, such as human body. The traditional 
near-field coils, e.g., lumped element L/C loop reso- 
nators and microstrip transmission line resonators 
(6,7), face challenges to provide high SNR gain in 
the deep areas of human body. To address this prob- 
lem, far-field coils or radiative antennas, e.g., dipole 
(8,9) and monopole (10), have been applied and 
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demonstrated stronger B, field and higher SNR gain 
over the conventional near-field coils in the areas 
deep inside human head and body (8,/0). 

When the monopole or dipole techniques are 
employed for densely-spaced arrays, e.g., 8-channel 
or 16-channle arrays for human head imaging, the 
distance of two adjacent coil elements is smaller 
than 10 cm (1/10 of the wavelength at 300 MHz). 
Therefore the close arrangement of coil elements 
might lead to non-ignorable electromagnetic (EM) 
coupling, which could decrease the SNR and paral- 
lel imaging ability. Several decoupling methods 
have been successfully applied for traditional loop 
and microstrip line coil, e.g., element overlapping 
(/1), transformers (/2), and interconnecting L/C net- 
works (/3). However, these methods face challenges 
and not readily feasible for monopole or dipole 
arrays due to their unique structures. To attain suffi- 
cient isolation among the radiative coil elements, 
the magnetic wall (MW) or induced current elimi- 
nation (ICE) decoupling technique (/4) has been 
proposed and successfully applied for monopole 
arrays (15) and dipole arrays (/6). 

In this study, we aim to theoretically analyze the 
MW decoupling for radiative coil arrays in ultra- 
high field MRI and derive the conditions required 
for achieving sufficient decoupling between radia- 
tive elements. Several theoretical methods, includ- 
ing the eigenvalue/eigenvector method, mutual 
impendence method, and odd/even mode method 
were employed in the analysis. Decoupling of a coil 
array is a process of modifying its impedance 
matrix and reducing all the off-diagonal elements 
to zero. For eigenvalue/eigenvector method, this 
can be regarded as a process of equalizing the 
eigenmode impedances. Mutual impedance method 
directly rearranges the impedance matrix to a diago- 
nal matrix. Compared with eigenvalue/eigenvector 
method, mutual impedance method is more straight- 
forward and easier to understand. However, it is 
limited by the complexity of the impedance matrix 
and is not suitable for analyzing arrays with high- 
count channels. Odd—even mode is another widely- 
used method for analyzing two-element arrays. This 
analysis assumes that any coupling between two 
coils can be decomposed into two basic modes, 
termed even and odd. In the even mode, the cur- 
rents on both coils have the same direction and 
amplitude. In the odd mode, the currents on the two 
coils have the same amplitude but opposite direc- 
tion. The coupling was totally reduced when the 
impedances in odd and even modes are the same. 
To validate the analytical results, transmission coef- 
ficient (S21), current distribution, and magnetic field 
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Figure 1 Illustration of two radiative elements and a 
decoupling element. V; and J; are voltage and current 
across the ports. Z; are the mutual impedances of the 
ports. Z; is the self impendence of the series circuit of 
the decoupling element. 


profile of MW decoupled radiative arrays were 
numerically studied and compared with those of 
radiative arrays without decoupling treatments. 


Il. THEORY 


Figure 1 shows a coil array consisting of two identi- 
cal radiative elements and a decoupling element 
which is placed symmetrically between the two 
radiative elements. The self and mutual impedances 
of the two radiative elements (Element 1 and Ele- 
ment 2) and the decoupling element (Element d) are 
represented by Z; and Z;;, respectively, where i Æj, 
{i,j} =1, 2, d. V; and J; are the voltage and current 
across the ports of the radiative elements and the 
decoupling element. The decoupling element 
includes a decoupling antenna and a series circuit. 
The self-impedances of the decoupling antenna and 
the series circuit are Zy and Zz, respectively. 

The current of the decoupling element is totally 
passive. Therefore, the voltage of the decoupling 
element is zero and the matrix notation of voltage 
and current is 


Zu Z2 Lal lh Vi 
Zr, Zo. Loa | | lo | =| V2 [1] 
Za Za Zaa | (Ta 0 
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Due to the reciprocity and the symmetry of Ele- 
ment | and Element 2, we can get 


Zai = Zia = Za Zza 
Zu =Zo2 


[2] 


Eigenvalue/Eigenvector Method 


In our previous study (/4), the two coil elements 
are decoupled when one of the eigenvalues of 
impedance matrix Z is equal to 0. It is worth noting 
that although the eigenvalue/eigenvector method is 
used for reactance matrix X, it could be used for the 
impedance matrix Z. The requirement of the decou- 
pling element from Ref. (14) is 
2 


Zz 
Zum [3] 
Z12 


Mutual Impendence Method 


From Eqs. [1] and [2], the matrix notation of the 
voltage and current of Element 1 and Element 2 
can simplified as (/7): 


Zi2— zia Z> zia 

Zdd Zdd 

The two coil elements were decoupled when the 

new mutual impendence matrix Z could be simpli- 

fied as a diagonal matrix. In other words, we 

require Zn- 4 =0. Thus the requirement of the 
decoupling element is 


Z 
Zua= £. [5] 


Odd/Even Mode Method 


For the odd mode, the voltages and currents of the 


two ports are: V, Vo=V3 h I, =I. Thus Eq. 
[1] can be written as: 
Zi Zi Lia I V 
Za Zn Zaļ||-I|=5|-V [6] 
Za Za Zaa] | la 0 


For the even mode, the voltages and currents of 
the two ports are: Vi = V2 = V; Jı = h} =I. Thus Eq. 
[1] can be written as: 


Zu Zi2 Zia I V 
Zn Z2 Ząaļ||ļ|IĪI|=]|V [7] 
Za Za Zaa] | La 0 


From Eqs. [6] and [7], the impedances of the 
odd-mode (referred to as Zə) and even-mode 
(referred to as Ze) can be derived as follows: 


V 
Z= T =Z -Z112 [8] 
2 
EE E E [9] 
I Zdd 


If Z,=Z,., the elements are decoupled. Thus, 
based on Eqs. [8] and [9], the decoupling condition 
can be derived as 


2 
Få 
Zaa= 14 [10] 


It is clear from Eqs. [3], [5], and [10] that the 
three theoretical methods gave the same require- 
ment of the decoupling element. Given that the 
decoupling element includes a decoupling antenna 
(Zaq) and a series circuit (Zz), the required imped- 
ance of Zz: 


2 
ZL= = —Lara [11] 
Z12 


Ill. SIMULATION 


A set of EM simulations using a full-wave EM solver 
(HFSS, ANSYS, Canonsburg, PA) were employed to 
validate the theoretical analysis described above (2). 
In this article, we only use monopoles as an example 
to demonstrate this concept. However, the basic prin- 
ciple will work for any coil, e.g., dipole array, since 
the method only relies on self and mutual impedan- 
ces. Two monopole elements with and without the 
decoupling element were modeled in simulation, as 
shown in Fig. 2(A,B). In the simulations, the length 
of the decoupling element was chosen the same as 
the monopole elements. Monopole elements were 
mounted on a cylindrical former with a diameter of 
25 cm. The angle of two monopole elements was 
varied from 25 to 50 degree with a step of 5 degree 
[Fig. 2(C)], which means the distance between the 
two monopoles was varied from 4.3 cm to 10.6 cm. 
The operating frequency is 297.2 MHz, which corre- 
sponds to the Larmor frequency of our 7 T MRI sys- 
tem. A cylindrical water phantom with an outer 
diameter of 16 cm and a length of 37 cm was placed 
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Figure 2 Three-dimensional (3-D) models of 2-channel monopole arrays without (A) and with (B) the decoupling ele- 
ment. C: Cross view of the two monopoles and the water phantom. Monopole elements were mounted on a cylindrical for- 
mer with a diameter of 25 cm. The angle of two monopole elements was varied from 25 to 50 degree with a step of 


5 degree. 


4.5 cm below the monopole elements. The EM 
parameters of the water phantom were set as: conduc- 
tivity ø = 0.59 S/m; relative permittivity ¢, = 78. The 
distance between the boundary and models was larger 
than A/2. Manual mesh was used to accelerate simu- 
lation convergence and the convergence condition AS 
was set to 0.002 to achieve more reliable results. 

As shown in Fig. 3, the procedure of the valida- 
tion is: 


1. The series circuit of the decoupling element is 
replaced with a port. Then we get self impedance 
of the decoupling monopole (Z 7) and the 
mutual impedances of monopole-monopole (Z12) 
and monopole-decoupling element (Zia). 

2. Secondly, the impedance of the serious circuit of 
decoupling element (Zz) is obtained by solving 
Eq. [11]. 

3. Finally, the equivalent lump component of Zz is 
employed in new simulation. S$, of the two 
monopoles, current distribution along the conduc- 
tors, and the magnetic field on the phantom were 
obtained and compared with the 2-channel array 
without decoupling treatments. 


IV. RESULTS 


Validation of S-Parameter 


Z12, Zig, and Zy,~ obtained from Step 1 of Fig. 3 
were listed in Table 1, from the second row to 
fourth row. The desired Zz; calculated from Eq. [11] 
was listed in the fifth column. The required resistan- 
ces of the decoupling elements (referred to as 
R(Z,)) are assumed to be neglectable because their 
values are very small. Thus only the imaginary 
parts of Zz (referred to as Im(Z;)) were replaced by 


lump components (capacitors or inductors), as listed 
in the sixth column. In order to validate the analysis 
described above, the equivalent lump components 
(sixth column of Table 1) were employed in the 
new simulation (Step 3 of Fig. 3). 

Figure 4 shows S>, of the two monopole ele- 
ments with and without decoupling elements. It can 
be observed that the coupling was reduced from 
about —8 dB to a sufficiently small value (<—35 
dB) when the predicated decoupling elements were 
employed. This indicates that the theory described 
above was reliable and can accurately estimate the 
circuit of the decoupling element. These results 
have also demonstrated that the assumption of 
neglectable R(Z,) is reasonable. Thus only reactance 
(capacitance or inductance) is needed and no Ohm 
loss is induced in the decoupling element. 


Validation of Current Distribution 
and H Field 


Figure 5(A) shows the current density distribution 
of two coupled monopoles when the left monopole 
was excited with | W while the right monopole was 
not excited. Due to the mutual coupling, the signal 
power was obviously coupled to the right monop- 
ole, especially when the angle equals to 20 or 25. 
By using the predicated decoupling elements, how- 
ever, the coupled signal became neglectable, as 
shown in Fig. 5(B). This comparison results were in 
agreement with the S>, results as described above. 
The current directions were indicated by black 
arrows in Fig. 5(A,B). The current of right monop- 
ole using MW decoupling was very small and thus 
its direction was not shown. 
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Setp3: Comparison of monopole arrays with (right) 
and without decoupling treatments (left) 
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Figure 3 The procedure of validating the theoretical analysis with a EM simulation study. 


Figure 5(C,D) shows the normalized H field on 
the water phantom in the transverse plane of 
monopole arrays without and with MW decoupling 
method. Clearly, the decoupling elements have 


shielding effect, making H field profiles more inde- 
pendent and weaker at the neighborhood of the two 
elements. As shown in previous studies (78), the shield- 
ing effect could benefit the MR imaging in terms of 


Table 1 Self and Mutual Impedances of the Monopole Elements and Decoupling Elements 


Angle (degree) Zı2 (Q) Zia (Q) Zra Z[11] (Q) Equivalent of Im (Zz) 
20 29.637 + j5.438 31.769 + j11.39 32.9 + j18.673 0.145 — 70.2816 1.9017 nF 
25 27.848 + j2.2095 31.087 + j9.3953 32.686 + j18.399 0.3036 — j0.0403 13.288 nF 
30 25.8 — j0.63357 30.271 + 77.5383 32.409 + 718.177 0.451 + 70.3192 0.171 nH 
33 23.551 — j3.0622 29.9362 + j5.8832 32.071 + j17.979 1.9913 + j1.4065 0.7536 nH 
40 21.215 — j5.0843 28.402 + j4.386 31.728 + j17.91 0.7114 + j1.6080 0.8615 nH 
45 18.86 — j6.7665 27.462 + j2.9963 31.507 + j18.248 0.7249 + j2.0418 1.094 nH 
50 16.437 — j8.0404 26.338 + j1.689 31.061 + j17.975 0.7168 + j2.9824 1.5979 nH 
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Figure 4 S2; plots vs. frequency of the two monopole elements with and without the magnetic wall decoupling method. 
The series circuits of the decoupling elements used in simulation were calculated from Eq. [11]. It is clear that the strong 
coupling could be reduced from —6 to —10 dB to a sufficiently small value (better than —35 dB). These results indicate 
that the theory described above was reliable and it can accurately estimate the requirement of the decoupling element. 


Figure 5 Surface current distributions of monopole arrays and magnetic field profiles on the water phantom. During the 
simulation, only the left port was excited with 1 W. A: Current distribution of the monopole array without decoupling 
treatments. B: Current distribution of the MW decoupled monopole array. By using the predicated decoupling element, the 
coupled current on the right monopoles can be well suppressed. Black arrows in (A) and (B) indicate the current directions 
along the conductors. C: Magnetic field profile of monopole array without decoupling treatments. D: Magnetic field profile 
of the MW decoupled monopole array. The magnetic fields near the right monopole become weaker when MW decoupling 
was employed, which is probably because the decoupling elements act as “magnetic wall” and show shielding effect. 
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SNR and parallel imaging performance. It can be under- 
stood that the shielding effect was generated by the 
induced current of the decoupling element. 


V. DISCUSSIONS AND CONCLUSION 


Magnetic wall/ICE decoupling method uses an inde- 
pendent resonator/antenna which has no physical 
connection to array elements to compensate for the 
current induced by EM coupling in array elements. 
In this study, this decoupling method applied for 
radiative coil array in ultrahigh field MRI was ana- 
lyzed by using the eigenvalue/eigenvector method, 
mutual impendence method, and odd/even mode 
method. The derived conclusions on decoupling 
conditions from all the three theoretical methods are 
consistent. These theoretical methods and the 
decoupling conditions were validated by a set of 
comparative EM simulations. 

The results of numerical analysis on transmission 
coefficient (S21) between array elements and the 
induced current distribution have demonstrated that 
the strong coupling (~—8 dB) could be reduced to a 
sufficiently small value (better than —35 dB) by using 
the decoupling condition obtained from the theoreti- 
cal analysis. The induced current of the decoupling 
element caused shielding effect and focused the mag- 
netic field of monopole element (Fig. 5). This focus- 
ing generates locally stronger and spatially diverse B, 
fields of antenna array elements, which might be 
advantageous for MR imaging in terms of SNR and 
parallel imaging performance. 

From Fig. 4, we can find that the decoupling per- 
formance still could be better than —30 dB even 
when the distance between monopoles is only about 
5 cm. This indicates that the ICE decoupling 
method is potential to design radiative coil arrays 
with high-count channels, e.g., 16 channels. Since 
the above derivation is purely based on self and 
mutual impedances, the decoupling element can be 
arbitrary and need not be of the same type. In other 
words, the decoupling element could be loop or 
microstrip line resonators as long as it meets the 
decoupling condition in Eq. [3]. 
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